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TRANSIENT  RESPONSE  OF  VOLUMETRIC  SOURCES 
AND  LINE  ARRAY  SYSTEMS 

First  viewgraph,  please  (1). 

' - - - \ 
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•  CONTINUOUS  LINE  ARRAY  SYSTEM 

•  FREQUENCY  FOCUSED  AND  SCANNED  ARRAY 

•  VOLUMETRIC  ARRAY 
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Viewgraph  1 

Although  considerable  work  has  been  devoted  to  the  analysis  of  the  steady 
state  response  of  line  arrays,  less  effort  has  been  devoted  to  a  systematic 
analysis  of  the  transient  response  of  these  arrays  for  wideband  signals.  The 
subject  of  this  talk  is  the  extension  of  the  impulse  response  technique  to 
free  field  discrete  line  array  systems  with  line  elements  and  time  delay 
beamformers.  We  will  examine  some  relatively  simple  discrete  and  continuous 
array  systems  before  showing  the  transient  response  of  a  complex  frequency 
focused  and  scanned  array  system  obtained  by  the  impulse  response  technique. 

We  will  then  show  how  the  technique  may  be  easily  extended  to  investigate  the 
field  of  a  volumetric  array  that  simulates  the  characteristics  of  a  transient 
parametric  array  system. 


—  Next  viewgraph,  please  (2). 
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Viewgraph  2 

Consider  the  problem  of  determining  the  time-dependent  pressure  at  a 
spatial  point,  r,  resulting  from  a  specified  spatial  and  temporal  velocity 
distribution  of  an  array  of  sources  in  an  ideal  fluid  of  density,  p,  and 
acoustic  propagation  speed,  c.  The  inhomogeneous  wave  equation  formulated  in 
terms  of  a  velocity  potential,  <*>,  and  source  term,  Q,  is  given  by  equation 
1.  The  pressure  may  then  be  obtained  from  the  velocity  potential  by  applying 
equation  2.  The  solution  to  the  wave  equation,  using  a  Green's  function 
approach  for  an  unbounded  medium  and  when  the  source  term  is  independent  of 
its  location,  is  given  by  equation  3. 

—  Next  viewgraph,  please  (3).  — 
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Viewgraph  3 

We  now  define  a  spatial  impulse  response  h  of  r  and  t.  This  response 
function  is  for  a  point  in  space  resulting  from  a  distributed  source  excited 
at  time  t0.  In  the  equation,  t0  is  an  excitation  time  and  may  be  related 
to  a  time-delay  beamformer.  The  quantity  R  over  C  is  related  to  the  problem 
geometry  and  is  simply  a  travel  time  between  the  source  and  spatial  point  of 
interest.  The  velocity  potential  may  then  be  expressed  as  a  convolution.  We 
can  then  write  an  expression  for  the  pressure  in  terms  of  a  convolution  of  a 
source  q  of  t  and  a  spatial  impulse  response  h  of  t. 

—  Next  viewgraph,  please  (4).  — 
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The  impulse  response  at  a  field  point  resulting  from  a  simple  source  is 
merely  an  impulse.  A  straightforward  extension  of  the  simple  source  to  a  line 
array  consisting  of  N  +  1  simple  sources  located  a  distance,  d,  apart  is 
illustrated  in  the  upper  part  of  this  viewgraph.  When  all  sources  are  excited 
simultaneously,  the  farfield  impulse  response  is  simply  a  series  of  impulse 
functions  beginning  at  a  time  corresponding  to  the  travel  time  to  the  nearest 
element.  The  train  of  impulses  ends  at  a  time  corresponding  to  the  farthest 
array  element  and  the  separation  between  each  impulse  is  simply  dependent  on 
the  angle,  e,  the  interelement  spacing,  d,  and  speed  of  propagation,  c. 

—  Next  viewgraph,  please  (5).  — 
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When  the  general  source  distribution  in  space  is  such  that  it  can  be 
represented  as  a  continuous  line  source,  a  closed  form  solution  for  the 
impulse  response  may  be  formed  by  a  linear  superposition  of  simple  sources. 

We  first  consider  the  case  when  the  field  point  is  in  the  farfield  and  outside 
the  region  defined  by  X  =  d,  as  shown  in  the  upper  section  of  viewgraph  5. 

The  impulse  response  in  this  case  is  simply  a  boxcar  function  beginning  at  a 
time  corresponding  to  the  nearest  end  of  the  source  and  ending  at  a  time 
corresponding  to  the  farthest  edge.  The  spatial  impulse  response  for  a 
complex  array  may  have  many  such  overlapping  boxcar  functions.  Because  it  is 
difficult  to  depict  overlapping  boxcar  functions,  we  have  decided  to  display 
the  time  derivative  of  the  impulse  response  when  appropriate.  A  sample  of 
this  is  shown  at  the  bottom  of  this  viewgraph.  A  positive  impulse  corresponds 
to  the  leading  edge  and  the  negative  impulse  corresponds  to  the  trailing  edge 
of  the  spatial  impulse  response  of  an  ideal  continuous  line  source. 

—  Next  viewgraph,  please  (6).  — 
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A 

TIME  DERIVATIVE  OF  IMPULSE 
RESPONSES  FOR  FIVE  ELEMENT  ARRAY 
WITH  AND  WITHOUT  BEAMFORMING 
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Viewgraph  6 


Consider  an  array  system  composed  of  line  elements  and  a  time  delay 
beamformer  as  shown  in  (a)  of  viewgraph  6.  The  individual  elements  are 
extended  sources  of  length,  L,  for  which  the  spatial  impulse  response  can  be 
represented  as  a  boxcar  function.  The  impulse  response  can  be  ideally  delayed 
by  the  beamformer.  As  a  simple  example,  consider  the  5-element  col  linear 
array  and  beamformer  illustrated  here  If  the  beamformer,  depicted  here  as 
parallelograms,  is  removed  by  setting  all  time  delays  (t-|  through  t5)  to 
zero,  the  derivative  of  the  transmit  impulse  response  at  the  field  point  is  as 
shown  in  (b).  The  strength  of  the  positive  and  negative  impulse  pairs  is 
determined  by  the  distance  between  the  center  of  the  element  and  the 
x-coordinate  of  the  field  position.  The  location  of  each  impulse  in  the  time 
domain  is  determined  by  the  element  location  and  spatial  extent,  L.  With  the 
field  point  located  directly  above  the  second  element,  the  impulses  from 
elements  1  and  3  are  collocated  and  simply  summed  using  linear  superposition 
as  indicated  by  the  stacked  arrows  in  (b).  When  the  beamformer  is  used  to 
focus  the  array  at  the  field  point,  the  temporal  spread  caused  by  the  element 
location  with  respect  to  the  field  point  is  removed.  The  temporal  spread 
caused  by  the  spatial  extent  of  the  elements  still  exists  as  in  (c).  That  is 
to  say,  as  expected,  the  directional  characteristics  of  the  elements  are  not 
removed.  Therefore,  if  the  total  time  spread  of  the  array  impulse  response 
(shown  in  (c))  is  small  compared  to  the  highest  frequency  component  of  the 
input  time  function  and  the  pulse  duration  is  greater  than  the  impulse  spread, 
a  steady  state  solution  approaching  an  ideal  array  output  will  be  obtained 
from  the  convolution  integral  or  impulse  response  approach. 
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—  Next  viewgraph,  please  (7).  - 
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CONCEPTUAL  DIAGRAM  OF  THIRTY-FIVE 
ELEMENT  COLINEAR  LINE  ARRAY 
AND  BEAMFORMER 
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Viewgraph  7 

The  transmit  transient  response  of  more  complex  col  linear  arrays  can  now 
be  examined.  Consider  the  ultrasonic  array  system  shown  here.  This  35-element 
array  has  an  interelement  spacing  of  3  mm  and  element  lengths  are  1.5  mm,  or 
one  wavelength  at  1  MHz.  The  array  beamformer  values,  tn,  are  set  such  that 
the  the  array  is  focused  at  the  point  indicated  by  the  circle  at  x  =  3  cm  and 
y  =  9  cm.  We  will  now  show  the  spatial  impulse  and  normalized  transient 
pressure  response  at  the  two  points  in  space  marked  by  asterisks. 

—  Next  viewgraph,  please  (8).  — 
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The  derivative  of  the  impulse  response  at  the  focal  point  is  shown  in  the 
upper  figure.  Note  that  the  beamformer  time  delays  have  reordered  the  impulses 
so  that  they  nearly  coincide.  The  remaining  temporal  dispersion  is  due  to  the 
finite  extent  of  the  radiators.  Large  differences  in  the  strength  of  the 
impulse  pairs  occur  because  the  effects  of  spreading  and  element  directivity 
are  included  in  the  impulse  responses.  The  transient  pressure  response  arising 
at  this  field  point  from  a  20  ys,  1  MHz  pulse  is  shown  in  the  second  part  of 
the  viewgraph.  We  see  that  the  excitation  waveform  is  recovered  and  a  steady 
state  solution  is  achieved.  The  spatial  impulse  response  and  corresponding 
pressure  response  due  to  a  similar  excitation  at  a  field  point  other  than  the 
focal  point  are  shown  in  the  lower  half  of  this  viewgraph.  We  see  the  spatial 
impulse  response  is  nearly  25  ys  in  duration  and  the  impulses  are  spaced 
unequally.  This  results  in  temporal  smearing  and  a  complex  interference  in 
the  pressure  response  function. 
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—  Next  viewgraph,  please  (9).  - 
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IMPULSE  RESPONSE  OF  CONTINUOUS  A 
LINE  SOURCE  AS  FUNCTION  OF  rj 
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Viewgraph  9 


For  simplicity,  the  previous  development  for  a  continuous  line  assumed 
the  source  was  uniformly  excited.  Now,  consider  a  case  of  the  more  general 
impulse  response  where  the  excitation  of  each  point  is  dependent  on  its 
position  and  the  speed  of  propagation,  v,  down  the  source.  This  translates 
into  a  delay  associated  with  the  excitation  function,  which  is  spatially 
dependent  on  the  source  location.  The  speed,  v,  could  be  associated  with  the 
speed  of  electromagnetic  propagation,  the  speed  of  a  propagating  wave  in  a 
different  medium,  or  a  well-designed  delay  line  for  a  beamformer. 


We  now  define  a  function,  n,  as  shown  in  the  upper  left-hand  corner  of 
viewgraph  9.  Let  us  examine  the  impulse  response  of  this  system  for  different 
values  of  n,  which  subdivides  the  field  into  regions  with  different 
characteristics  for  the  impulse  responses.  When  n  >  0,  the  leading  edge  of 
the  impulse  arises  from  the  source  point  farthest  from  the  field  point  and  the 
trailing  edge  location  is  dependent  on  the  values  of  v,  c,  and  0.  Thus,  the 
path  directly  through  the  medium  is  the  fastest.  When  n  <  0,  the  reverse  is 
true  and  the  leading  edge  is  dependent  on  v,  c,  and  e,  so  it  is  faster  to 
travel  down  the  source  than  through  the  medium  to  the  field  point.  When 
n  =  0,  then  sin  n  =  c/v,  and,  along  this  radial  in  the  farfield,  all 
elemental  point  sources  contribute  at  the  same  time  forming  a  single  impulse 
at  t  =  R/C.  The  significance  of  the  result  is  that  the  proper  delays  required 
for  beamforming  can  be  obtained  directly  from  the  array  impulse  response 
function.  We  now  examine  the  pressure  responses  for  this  system. 


—  Next  viewgraph,  please  (10).  - 
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Viewgraph  10 

Consider  an  example  where  the  array  length  is  0.2  m  and  the  delay  line 
speed  is  2000  m/s  (as  in  (a)).  The  three  field  points  of  interest  lie  on  the 
arc  of  a  10  m  circle  at  bearings  such  that  the  three  possible  ranges  of  n 
are  examined.  The  transient  pressure  response  resulting  from  a  20  ys  CW 
excitation  pulse  at  each  of  the  three  field  points  is  presented  in  viewgraph 
10.  When  n>0,  the  time  (tfa)  at  which  a  disturbance  first  reaches  the 
field  point  is  associated  with  the  end  of  the  array  located  at  the  coordinate 
system  origin.  The  transient  duration  is  about  55  ys  and  a  steady  state 
pressure  response  is  never  achieved  In  (b),  the  pressure  function  appears  to 
be  a  near  replica  of  the  input  signal  as  expected  when  n  =  0.  There  is  a 
transient  effect  near  the  beginning  and  end  of  the  response,  but,  for  the  most 
part,  this  figure  represents  a  steady  state  solution.  The  arrival  time  of  the 
initial  disturbance  is  as  expected,  i.e.,  the  same  as  in  the  n  >  0  case. 
Furthermore,  it  should  be  noted  that  the  ratio  of  the  pressure  amplitudes 
shown  in  (a)  and  (b)  is  nearly  60:1.  When  n  <  0,  the  initial  disturbance 
arises  from  the  end  of  the  array  nearest  the  field  point.  This  means  it  is 
faster  for  a  disturbance  to  travel  through  the  array  before  entering  the 
medium  than  to  travel  solely  in  the  medium.  The  transient  response  begins  at 
6652  ys  and  lasts  about  35  ys  so  that  for  the  n  <  0  case  the  end  of  the 
transient  is  the  same  as  the  beginning  of  the  transient  when  n  >  0.  As  in 
the  case  for  n  >  0,  the  amplitude  of  the  transient  is  substantially  less 
than  in  the  case  when  n  =  0. 
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—  Next  viewgraph,  please  (11).  - 
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CONCEPTUAL  DIAGRAM  OF  FREQUENCY 
FOCUSED  AND  SCANNED  LINE  ARRAY 
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Viewgraph  11 

A  further  application  of  the  impulse  response  technique  to  a  more  complex 
array  system  is  worthwhile.  Consider  an  electronically  focused  and  scanned 
line  array.  This  type  of  array  system,  used  in  acoustic  imaging,  has  been 
described  by  Kino  and  coworkers.  They  have  shown  a  time-dependent  frequency 
signal  sent  through  a  delay  line  can  be  used  for  the  purpose  of  electronically 
focusing  and  scanning.  This  viewgraph  illustrates  the  concept  in  an  array 
system  of  this  type.  The  initial  excitation  time  of  each  element  is  delayed 
an  amount  determined  by  the  delay  line  speed  and  element  position.  The 
element  remains  energized  for  the  pulse  length.  Normally,  the  pulse  length  is 
less  than  the  time  required  for  a  disturbance  to  transit  the  delay  line.  In 
effect  this  results  in  a  subset  of  elements  energized  at  any  instant  of  time. 
This  subset  forms  an  active  array  that  transits  the  array  elements  at  speed, 
v,  and  is  the  basis  of  a  scanned  array.  For  focusing,  a  linear  frequency 
modulation  (LFM)  chirp  signal  is  used  in  practice  as  it  is  easily  generated. 
The  selection  of  the  center  frequency,  which  is  dependent  on  the  interelement 
spacing.  Ax,  and  delay  line  speed,  v,  determines  the  azimuthal  angle  of  the 
maximum  beam  while  the  LFM  slide  constant  determines  the  focal  range.  The 
impulse  response  technique  can  be  readily  used  to  analyze  the  behavior  of  such 
complex  array  systems. 


—  Next  viewgraph,  please  (12). 
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FARFIELD  SPATIAL  IMPULSE  AND 
TRANSIENT  PRESSURE  RESPONSE 
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Viewgraph  12 
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Consider  now  an  array  of  100  equispaced  elements  when  the  focal  range  is 
large  with  respect  to  the  entire  array  length,  L.  The  impulse  response  for  a 
typical  field  point  is  shown  at  the  top  of  the  viewgraph.  For  this  case,  the 
amplitude  of  each  impulse  will  be  nearly  constant.  Furthermore,  the  arrival 
time  difference  between  each  impulse  is  nearly,  but  not  exactly,  the  same. 
Thus,  for  a  sample  point  in  the  farfield,  the  focusing  process  is  achieved 
with  a  CW  pulse  traveling  down  the  array.  The  instantaneous  transient 
pressure  response  obtained  using  the  impulse  response  technique  and  the 
indicated  variables  are  shown  in  the  lower  half  of  the  viewgraph.  The  sine 
dependence  on  time  is  expected.  It  can  be  noted  that  the  ratio  of  mainlobe  to 
first  sidelobe  is  very  near  the  expected  value  of  4.7,  or  equivalently  13.5 
dB.  Also  of  interest  is  the  slight  asymmetry  of  the  response  in  time  and 
amplitude  because  of  aberration  induced  by  the  pulse  traveling  down  the  delay 
line. 


—  Next  viewgraph,  please  (13).  — 
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A 

SPATIAL  IMPULSE  AND  TRANSIENT 
RESPONSE  OF  FREQUENCY  FOCUSED 
AND  SCANNED  LINE  ARRAY 
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Consider  now  a  case  when  the  field  point  is  relatively  close  to  the 
array,  but  still  at  a  range  such  that  each  element  behaves  as  a  simple  source. 
The  strength  of  the  impulse  response  arising  from  each  element  is  quite 
simple,  but  its  location  in  time  relative  to  other  elements  is  more  complex. 
The  impulse  response  for  a  100-element  array  at  a  focal  range  of  10  cm  is 
depicted  in  the  upper  part  the  viewgraph.  In  this  particular  case,  the  field 
point  is  located  over  the  array  midpoint.  Thus,  though  difficult  to  see, 
there  are  as  many  impulses  to  the  left  as  to  the  right  of  the  maximum  value. 
However,  the  temporal  spread  about  the  maximum  value  is  not  symmetrical.  The 
exact  frequency  required  to  obtain  the  constant  phase  values  needed  for 
focusing  can  be  readily  obtained  by  determining  the  time  between  impulses. 

For  a  given  pulse  length,  only  a  portion  of  the  array  is  excited  at  any  given 
time.  If  the  exact  frequency  content  required  for  focusing  is  replaced  by  an 
approximation  (as  in  the  case  with  an  LFM  chirp),  then  a  reasonable  facsimile 
of  ideal  focusing  can  be  achieved.  When  the  LFM  signal  dictated  by  this 
technique  and  containing  a  center  frequency  of  2  MHz  is  convolved  with  the 
impulses  (as  shown  in  the  upper  figure)  and  differentiated,  the  pressure  field 
shown  in  the  lower  figure  is  obtained.  It  is  clear  that  the  array  system  is 
focused  at  the  spatial  point  of  interest  about  32  ys  after  the  initial 
arrival . 

There  are  several  other  features  in  the  transient  pressure  response  that 
can  be  explained  on  the  basis  of  the  impulse  response  function.  Note  from  the 
impulse  response  in  the  upper  figure  that  the  temporal  samples  are 
nonuniformly  spaced  in  time.  This  is  analogous  to  unequal  element  spacing.  It 
is  well  known  from  antenna  design  theory  that  the  radiation  pattern  of  an 
aperiodic  spatial  array  may  be  described  in  three  parts.  The  maximum  response 
region  exhibits  characteristics  for  which  the  array  is  designed.  This  is 
followed  by  a  clean  sweep  region  and  then  a  region  of  moderately  high  levels 
called  the  plateau  region.  These  regions  are  then  analogous  to  the  time 
domain  response  observed  in  the  lower  figure  with  a  few  complications.  Observe 
that  there  is  a  maximum  response  at  approximately  32  *is;  there  are  regions 
of  very  low  amplitude  adjacent  to  the  maximum  response;  and  there  is  the 
resemblance  of  plateau  regions  near  the  beginning  and  end  of  the  transient 
response.  The  picture  is  somewhat  complicated  because  in  this  instance  the 
amplitude  response  is  dependent  on  time,  and  a  moving  pulse  leads  to 
nonsymmetrical  results  because  of  aberration. 

—  Next  viewgraph,  please  (14).  — 
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TRANSIENT  RESPONSE  OF  VOLUMETRIC 
ARRAY  USING  IDEALIZED  VOLUME 
SOURCE  DISTRIBUTIONS 
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Viewgraph  14 

As  a  final  case  of  interest  to  illustrate  the  usefulness  of  the 
generalized  impulse  response  approach,  we  note  that  the  technique  is  readily 
applicable  to  evaluate  the  characteristics  of  the  acoustic  transient  fields 
that  are  generated  by  idealized  secondary  source  distributions,  which  are 
broadband.  A  specific  example  of  interest  is  the  case  of  a  projector  or 
transmitting  array  that  is  generating  a  highly  directive  finite  amplitude 
broadband  signal  propagating  along  its  axis  as  illustrated  in  the  viewgraph. 
A  theory  to  evaluate  the  on-axis  transient  pressure  from  such  a  transient 
parametric  array  was  originally  presented  by  Berktay  and  was  later 
experimentally  verified  by  Moffett  and  co-workers. 

—  Next  viewgraph,  please  (15).  — 
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Viewgraph  15 

Without  presenting  the  details  of  the  development,  which  are  omitted  here 
because  of  time  constraints,  it  is  noted  that  the  characteristics  of  transient 
pressure  fields  resulting  from  transient  parametric  arrays  can  be  simply 
evaluated  via  the  present  impulse  response  approach.  The  edge  diffraction 
effects  are  ignored  in  the  source  term,  q,  which  is  represented  as  a  spatially 
decaying,  one-dimensional,  propagating  wave  in  the  cylindrical  volume  defined 
by  the  aperture  as  its  base.  The  term  W(X3)  in  the  source  term  denotes  the 
weighting  of  the  source  distribution  in  the  X3  direction,  which  for  a 
transient  parametric  array  would  be  exponential  in  nature.  The  farfield 
time-dependent  pressure  at  any  point  resulting  from  the  idealized  source 
distribution  can  be  expressed  as  a  convolution  of  the  source  function  s(t) 
with  the  impulse  response  hv(r,t).  Finally,  in  the  farfield  it  is  noted 
that  the  volume  impulse  response  can  itself  be  expressed  as  a  convolution  of 
two  impulse  response  functions  as  denoted  on  the  viewgraph.  These  latter 
impulse  responses  correspond  to  the  impulse  responses  for  an  equivalent  line 
source,  neglecting  the  aperture  effect  and  for  an  aperture  impulse  response, 
neglecting  the  line  effect.  It  is  noted  that  the  convolution  of  the  impulse 
responses  to  obtain  hv(r,t)  is  the  time  domain  equivalent  of  the  usual 
product  theorem  in  array  analysis  and  could  have  been  anticipated  apriori. 

—  Last  viewgraph,  please  (16).  — 
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•  THE  IMPULSE  RESPONSE  TECHNIQUE  IS  A  VIABLE  APPROACH 
FOR  ANALYZING  COMPLEX  ARRAY  BEAMFORMING  SYSTEMS 
SUBJECT  TO  STEADY  STATE  AND  TRANSIENT  EXCITATIONS. 

•  THE  TIME  HISTORY  OF  SIGNALS  GENERATED  OR  RECEIVED 
BY  SUCH  ARRAY  SYSTEMS  CAN  BE  STRONGLY  AFFECTED  BY 
THE  SYSTEM. 


•  THE  IMPULSE  RESPONSE  TECHNIQUE  CAN  BE  USED  TO 
EVALUATE  TRANSIENT  RADIATED  FIELDS  FROM  IDEALIZED 
VOLUMETRIC  SOURCE  DISTRIBUTIONS. 
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Viewgraph  16 

An  approach  has  been  developed  to  analyze  the  transmit  response  of 
complex  arrays  with  time  delay  beamformers.  This  approach  is  based  on  an 
extension  of  the  spatial  impulse  response  developed  by  Stepanishen. 

The  impulse  response  technique  is  a  viable  approach  to  the  formidable 
problem  of  analyzing  array  systems  subject  to  transient  excitations.  The 
technique  yields  a  spatial  impulse  response  that  is  dependent  on  the  geometry 
of  the  problem  and  beamformer.  The  impulse  response  can  then  be  convolved 
with  a  variety  of  signal  waveforms  to  perform  an  array  systems  analysis  as  a 
function  of  waveform  design.  This  approach,  therefore,  can  be  used  to  conduct 
a  rather  complete  systems  analysis. 
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